MODAL ANALYSIS OF SOME STRUCTURES USING 

CIRCLE - FIT METHOD 


by 

H. R. NIRANJANA MURTHY 



MU ^ 



V\0^ 


DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 


FEBRUARY, 1988 



MODAL ANALYSIS OF SOME STRUCTURES USING 

CIRCLE - FIT METHOD 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


by 

H. R. NIRANJANA MURTHY 


to the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

FEBRUARY. 198S 



13 

CtNTR 

I i 

Jcc. No. 


APR 1389 

'■L 

T , . A • 

Aloiliz 


lY) £ ■- f1(S(S- M - M (7^ - McO 



CERTIFICATE 



This is to certify that the thesis entitled 
'modal analysis of some STRUCTURES USING CIRCLE-FIT 
method' by H.R. Niranjana Murthy is a bonafide 
record of work done by him under ray guidance and 
supervision for the award of the degree of Master of 
Technology at the Indian Institute of Technology, 
Kanpur, The work carried out in this thesis has not 
been submitted elsewhere for the award of a degree. 


( Dr. BHUPINDER PAL SINGH ) 

Assistant Professor 
Department of Mechanical Engineering 
I.I.T. Kanpur -20 801 6 



ACKNOa'LEDGEMENT 


11 


I acknowledge with sincerity and gratitude, 
the guidance provided by Dr. Bhupinder Pal Singh in 
sorting out various problems through out the course of 
this thesis. I thank him for his contributions made by 
way of timely advices and criticisms . 

I would like to express ray sincere gratitude 
and thanks to Dr. V. Sundararajan for suggesting the 
topic of the thesis. 

I thank Mr. H.V.C. Srivastava for his careful and 
neat typing of this thesis. 


H.R, Niranjana Murthy 


Feb. 25, 1988 
I . I .T . Kanpur 



iii 


CONTENTS 

Page 

CERTIFICATE i 

ACKNOWLEDGEMENT ' 

CONTENTS iii 

LIST OF TABLES i^ 

LIST OF FIGURES ^i 

NOMENCLATURE vii 

ABSTRACT ' ix 

CHAPTER-I INTRODUCTION ^ 

1 .1 REVIEW OF PREVIOUS WORK 2 

1 .2 OBJECTIVES AND SCOPE OF PRESENT WORK 3 

CHAPTER-II CIRCLE-FIT METHOD ^ 

2.1 SODF STRUCTURE ^ 

2.2 MDOF STRUCTURES ^0 

2.2.1 Free Undamped Vibrations 10 

2.2.2 Free Damped Vibrations 12 

2.2.3 Forced Damped Vibrations 13 

2.3 CIRCLE -FIT METHOD 15 

CHAPTER-III RESULTS AND DISCUSSION 

3.1 MODELLING METHODOLOGY ^4 

3.2 RESULTS AND DISCUSSIONS 


REFERENCES 


53 




LIST OF TABLES 

iv 

Table No. 

Title 

Page 

3.1 

Comparision for a discrete system 

(2x2) (Structural damping) 

29 

3.2 

Comparision for an axial rod (2x2) 

(structural damping) 

31 

3.3 

Comparision for an axial rod (3x3) 

(structural damping) 

32 

3.4 

Comparision for a beam (4x4) 

(structural damping) 

34 

3.5 

Comparision for a beam (4x4) 

(Structural damping) 

35 

3.6 

Comparision for a beam (6x6) 

(Structural damping, b =0.01) 

36-37 

3.7 

Comparision for a beam (6x6) 

(Structural damping, b = 0.05) 

38-39 

3.8 

Comparision for L -Frame (6x6) 

(Structural damping, b=0.01) 

41 -42 

3.9 

Comparision for L-Frame (6x6) 

(Structural damping, b=o.05) 

43-44 

3.10 

Comparision for F -Frame (12x12) 

(Structural damping) 

45-46 



Table No. Title Page 

3.11 Comparision for a discrete system 
(2x2) (Viscous damping) 

3.12 Comparision for an axial rod (2x2) ^9 

(Viscous damping) 

3.13 Comparision for a beam (4x4) 50 

(Viscous damping) 

3.14 Comparision for L-Frame (6x6) 51-52 

(Viscous damping) 



Vi 


LIST OF FIGURES 


Figure No , 

Title 

Page 

2.1 

SODF Structure 

5 

2.2 

Nyquist plot for SDOF system 

7 

2.3 

Nyquist plot for MDOF structure 

17 

2.4 

Properties of modal circle 

19 

2.5 

Plots at w = w 

r 

22 

3.1 

Two DOF discrete system 

28 

3.2 

Axial rod, finite elements 

30 

3.3 

Beam, finite elements 

33 

3.4 

L- Frame structure 

40 

3.5 

F- Frame Structure 

45 



vii 


[C ] 

r^jk 

f(t) 


[H ] 


[M ] 


NOMENCLATURE 

Viscous damping coefficient 
Viscous damping matrix 
Diameter of modal circle 

Excitation force 

Hysteretic damping coefficient 

Hysteretic damping matrix 

Imaginary part 

Coordinates 

Stiffness 

Modal or generalized stiffness 

Stiffness matrix 

Mass 

Modal or generalized mass 
Mass matrix 

Numter of DOF in MDOF system 
Mode number 

Real part 

Undamped system natural frequency of SDOF 
Undamped system natural frequency 



w 


Excitation frequency 
x Displacement 

[Y((j)],y^j^ Mobility matrix, function 

. Velocity (first time derivative) 

.. Acceleration (second time derivative) 


[ a(w Receptance matrix, function 

P ,3' Constants 


r , 7 ' 


Constants 


T) 

r 




l.*]r 

[f] 



P 


t,Q 


Hysteretic damping loss factor 
Viscous damping ratio 
Eigen vector /mode shape 
Modal matrix 

Mass -normalized eigen vector /mode shape 

Mass -normalized modal matrix 
Square matrix 
Column matrix 

Row matrix 

Principal coordinates 
Angles 



lx 


ABSTRACT 


Modal analysis of some discrete and continuous 
structures has been performed in the present work. Circle fit 
method is adopted for the analysis. Analysis has been done for 
structures with viscous /structural proportional damping. 

The results of the identified system parameters (mass and 
stiffness distributions) have been compared with the 
actual values. The close agreement between identified and 
actual values is shown. 



CHAPTER -I 


I PRODUCTION 


The dynamic behaviour of real life structures are 
difficult to be analysed theoretically because of their 
complex geometry and uneven distribution of spatial parameters 
(mass, damping, stiffness). Even for the cases for which 
analytical and numerical solutions are possible, the nature of 
simplifications made necessitates an experimental verification. 
Modal analysis is one such method which is increasingly being 
employed for identification and validation of dynamic 
properties of physical systems. 

Modal analysis involves the testing of structures with 
the aim of obtaining a mathematical description of their 
dynamic and vibration behaviour. The mathematical description 
could be an estimation of natural frequency and damping factor 
in some cases and the full evaluation of mass, stiffness and 
damping distribution in other cases . The method involves 
test -structure excitation to obtain the modal responses like 
receptance, mobility, etc., and the extraction of system model 
from the test data. The entire procedure is sometimes also 
referred to as modal testing. 



1 .1 REVIEW OF PREVIOUS WORK 


A good review of the literature on modal analysis is 
given by Kundra [1 ] . Vepa [2] has recently reviewed modal 
identification methods and presented some of the inadequacies 
of the method. The topic of modal testing is neatly dealt with 
in a text by Ewins [3], 

One of the earliest and pioneer paper in the development 
of the subject is of Kennedy and Pancu [4 ] . They presented a 
circle fit method for the system identification. Bishop and 
Gladwell [5] presented a state of the art theory of resonance 
testing which at that time was considerably in advance of its 
practical implementation. Salter [6] presented a relatively 
non-analytical approach to the interpretation of the measured 
data. Salter’s approach, even though involving a good skill on 
the part of exp erimeniEr, gave a considerable insight into the 
vibration of structures . 

After 1970, owing to the great advancement in instrumen- 
tation, considerable progress has been achieved. A bibliography 
of several such works is given in [7,8 ] . A few of the important 
works are mentioned here. Flannelly et al [ 9 ] presented a theory 
to extract a complete physical model from multipoint sinusoidal 
excitation data. They assume that the structure has as many 
degree -of -freedom as the number of excitation points. Kundra [1 ] 
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has employed "this method in his work. He presented the theory 
for the lumped mass matrices. Berman [10 ] , in his paper, 
proposed a procedure which made use of the data obtained by 
the single point excitation. This approach makes the use of one 
column of mobility matrices at different frequencies, 

1 ,2 OBJECTIVE AND SCOPE OF PRESENT WORK 

The aim of the present work is to identify the modal 
parameters for discrete as well as continuous structures with 
viscous /structural proportional damping by modal analysis. Mass 
matrices for continuous structures are identified as consistent 
mass matrices. 

In the second chapter the theoretical aspects of the 
present method are discussed. The circle fit method, which is 
employed in the present work, is explained in detail. 

In the third chapter, the identified system parameters 
like mass and stiffness distributions are compared with the 
actual values both for discrete and continuous structures with 
different dampings mentioned earlier. Conclusions are presented 
at the end . 
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CHAPTER- II 

CIRCLE -FIT METHOD 

This chapter describes the circle fit method of modal 
analysis for multi degree of freedom (MDOF) structures. To 
understand the theory of this method for MDOF structures, it is 
worth to review the basics of vibration of single degree of 
freedom (SDOF) structure and MDOF structures, 

2.1 SDOF STRUCTURE 

The SDOF structure is shown in Fig. 2.1 , It consists of 
mass m, spring of stiffness k and massless damper with coefficient 
of viscous damping c or coefficient of hysteretic damping h, 
f(t) is the excitation force applied to the mass and x is the 
response of this mass , 

Its equation of motion with viscous damping is 
mx + cx + kx = f(t) =Fe^^^ .... (2.1) 

Assuming its solution as x^Xe^*^^, one gets the response as 

X = — I .... (2.2) 

(k-m m)+3wc 

Thus receptance for SDOF structure is 


Or(to) = ^ - 

1 

.... (2.3) 

X 

(k-Ki) si)+oo)C 


and its phase is 



/a(a)) = tan**^ 

o 

( -coc/k-oj'^m) ) 

(2.4) 


And its mobility is 
Y(m) = f = 


ju a((i!) 


(2.5) 
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It may be noted that receptance and mobility are complex 
quantities. Viscous damping model does not closely represent 
the real structures. Most of these structures when subjected to 
cyclic loading exhibit a stress -strain relationship which is a 
closed loop known as hysteresis loop. The energy dissipated per 
cycle due to this internal friction is proportional to the area 
of the hysteresis loop. Thus this damping is known as hysteretic 
or structural damping. This internal friction (structural 
damping) is found to be independent of frequency and is propor- 
tional to the displacement over a significant frequency range 
[11 "]. Since energy is dissipated, this damping force must 
be in phase with the velocity. Thus for simple harmonic motion, 
the damping force is given by jhx, where h is hysteretic damping 
coefficient. Equation of motion (2.1) for this hysteretic 
damping becomes 

mx + jhx + kx = Fe^^"^ .... (2.6) 

Receptance for this case is 

a(w) = I = 1 (2.7) 

(k-co m+ jh) 

For modal analysis purposes the receptances and mobilities, 
which are complex quantities arepref erably given in the Nyquist 
or Argand plots. Mobility response for viscously damped case 
and receptance response for hysterically damped case are given 
in Fig. 2.2 and are exact circles, proved as follows: 



(b) Receptance response with structural damping 


Fig. 2.2 Nyquist plot for SDOF system. 
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Equation (2.5) gives the mobility response as 


Y(<^) 


Jid. 


(k-Ci) m)+ juc 


.... ( 2 . 8 ) 


This can be rearranged as 
Y(a)) = 




'S o o' d 

(k-co m)'^+(coc)'^ 


a3(k-co m) 
(k-cci^m) ( CDC ) ^ 


or 


Y(w) = Re + j Im 


(2.9) 


It can be seen that 

(Re “ ■^)^+ Im^= (■^)^ .... (2.10) 

Thus the mobility response plot for viscously damped SDOF 
structure is an exact circle with radius (1/2c) and centre as 
(Re= Im= O) , Fig. 2.2(a) . 

Equation (2.7) gives the receptance response as 
a(w) = J 

k-u m + jh 


This is rearranged as 
a(w) 


(k-(o^m) 


_ = K A + j 

(k-to'^m) + h 


(-h)_ 

(k-w^m)^ + h^ 


or 


a(w)= Re + j Im 


.. (2.11) 


Real and imaginary parts are related as 

^)2= (^)2 


(2.12) 
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Thus the receptance response plot for structurally damped 
SDOF structure is also an exact circle with radius as (1 /2h) 
and centre as(Re=0,Im= - * Fig. 2.2(b). It may be noted that 

receptance and accelerance Argand plots for viscously damped 
case, and mobility and accelerance Argand plots for hysteretically 
damped case, are not exact circles but are circles like. 


The next step is to find the expression for damping by 
analysizing the modal circles. Considering the viscously damped 
case, it is observed from the Fig. 2.2(a) and the Eqn. (2.9) 
that. 


tanr = tan (■^) 


k-m^m 

me 


-2 2 
m -m 


(2.13) 


If one takes two points 90°, below natural frequency 

and Q= 90° above natural frequency as in Fig. 2.4(a), one can 

3 . 

write the following using Eqn. (2.13) 


1 = 


1 = 


2 -2 
CO -m^- 

a 

cwT/nT 


Adding above equations one gets 






m 


(2.14) 


Extending the same procedure to structurally damped 
system, one can write the following, by considering Eqn. (2.11) 
and Fig. 2.2(b) 



tanx 


h/m 

—2 2 

and tan 0/2 = tan(90-T) = 


(2.15) 


If one takes two points, one ©|^ = 90° below natural frequency 
Q and other ©^^=90° as shown in Fig. 2.4(b) above natural 
frequency, then using Eqn. (2.15), one obtains 



h'/m 



1 


Adding above equations, one gets 


2 , , 2 ^h ohw 

“a " “b - 


or 


h 

k 


4 - 




2w^ 


( 2 . 16 ) 


2.2 MDOF STRUCTURES 

In this section basics of vibration of MDOF are reviewed. 
2.2.1 FREE UNDAMPED VIBRATIONS 

The equations of motion for an undamped case, expressed 

in matrix form are 

# 

[m] {x} + [k] [x} = to} (2.17) 

Here brackets [ ] stand for square matrices and braces } 
stand for column matrices. Assuming [x} = (x) e^^^, the above 

equation reduces to 

[K] [X] = [M] (X] 


(2.18) 
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which is an eigen equation. 

Solving above equation, one gets eigen values (5^) and 
eigen vector matrix [f ]. 

The eigen vectors possess a very important and useful 
property known as Orthogonality property [12]. This is not 
an ordinary orthogonality, but an orthogonality with respect to 
mass matrix [M] as well as to stiffness matrix [K], and is 

Js [m] = 0 V 4 S (2.19) 

Js [K] Wr =0 r ^ s 


Here brackets |_ J stand for row matrix. 

The above equations show 

p f [H ] [? ] = 

[? f [K 1 [5 ] = 

where and are diagonal matrices . 


( 2 . 20 ) 


Eigenvectors [ ^ ] can be normalized as [^ ] , 

[Oi] =[ f - 7 .... (2.21) 


Then 




[f] [K] [0]- rsjq 


... ( 2 . 22 ) 
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2.2,2 FREE DAMPED VIBRATIONS 

All structures have some damping. For analysis purpose 
it is modelled as viscous, hysteretic, coulomb, aerodynamic etc. 
Generally it is difficult to determine which model is more 
appropriate for a given structure. Analysis becomes quite simple 
if viscous and hysteretic damping are proportional to stiffness 
and mass, i.e. 


[C] = ]B'[k] + -y'CM] 
[H] = P [K] + r [M] 


.... (2.23) 


This assumption has been used in this work. 


Equation of motion for viscous damping becomes 


[M] (x) + (^'[k] +‘r'[M] ) (x] + [k]{x] = [o] .... (2.24) 

Using coordinate transformation [x] = [t] (p) .... (2.25) 

and premultiplying by [T the Eqn, (2.24) becomes 

{ P] + ) (pj (p) ={o] 

( 2 . 26 ) 

The r-th equation is 

®r P'r ^r^ P^ + Pr = ^ 

_ . _2 

r T^r 


or 


0 


(2.27) 
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Equation of motion for hysteretic damping is 

[m ] (x) +oO [k] + y[n ]) [x] + [k ] {x} = [o] 

(2.29) 

\ 

Proceeding as above, this equation reduces to 

(1+j b^)Pj,=0 .... (2.50) 

where ,h = ^+r/C^ .... (2.31) 

r m^ r r 


2.2.3 FORCED DAMPED VIBRATIONS 

Equation of motion, for viscously damped case becomes 
[m] (x] +0'[k] +t'[m] ) (x) + [k] (x) = (F} 

.... (2.32) 

Assuming [x} = (X] one gets 

( [k] -0)^ [m] +0 w(3' [k] +7'[m] ) ) (X) = (Fj .... (2.33) 

or (XJ =[ a(m)] [Fj 

where 

[a(w) ]“^ = ([ k] -0)^ [m] + ow(/3'[k] +0^ '[m] )) 

(2.34) 

For modal analysis, it is better to change this [a(o)^ 
matrix in terms of modal properties. Premultiplying [a(co) ]” 
by [ff and post-multiplying by [ f ] to obtain 
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[^] [a(w) ]i [? ] = ( ['■m£r~-]+ jw( 3 ' L^k^l + T" ) ) 

Taking inverse on both sides, one gets 

,-1 


[ T ] [ a(w)] [ T ] ^ =( ["k^] [■^ni^>j(oO '[^kp^]+7 "["m^])) 

_ rn 

Again premultiplying by [ ^ ] and post^nultiplying by [ ^] gives 

[ a(o.) ] = [f ] ( [-k^] V '[-k^l +7 ) ) "'' [t ] 


T 


Typical coefficient of this matrix is 


.... (2.35) 


N 

a., (w) = L 


dk 


'i' . W 

-r->L.r-k. 


r=1 (k^-o)^m^) + OCo(g'k^+r'in^) 


or 


N 


a., (w) = E. 
jk^ ^ r=1 


^ ■ 0, 

r-J.,..£.!jL.... 

(c3^-o)^)+jco(3'co^+7' ) 


(2.36) 


where 0. is the j-th element of r-th mode . 

r j 

Equation of motion for hysteretic damping is 

[M](x}+a(3 [k] +7 [M]) (x}+[k1(x}= (f} (2.37) 

-1 


Proceeding as above [ a(^) ]*" matrix in this case is 
[a(w)]"''= ([K ] [m] +j(3 [K ]+7 [m])) 


(2.38a ) 


and its typical coefficient is 
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N 




r^.i r^k 


r=1 mj,)+d(3k^+'3'm^) 


or 


N 


d . o{. 

r^n r^k 


Ct -; V ( ^ ) “ ^ p p p 

r=1 (5^-o)^)+jOw>'y) 


(2.38b ) 


2.3 CIRCLE-FIT MEl'KOD 


This section describes the circle fit method for 
extracting natural frequencies, damping coefficients, mass and 
stiffness coefficients both for viscously damped and structurally 
damped cases . Firstly the effect of off -resonant contribution 
is briefed for structural damping case (same is true for viscous 
damping) , Secondaly the procedure for modal parameter extraction 
is discussed. 

Response of MDOF structures is given by Eqn,(2.33) 

{X} = [a(o)] (f) 


where 


ajk(“) = 


N 

E 

r=1 


r^.i r’^k 


(Cj-03^)+j(j3Q^+T) 


.... (2.38) 


Equation (2,38) can be expressed as 


r^k 

a-i (w)- — 5 — ^ ^ 5 




N 

s=1 (Cj-<»)^)+i( 35^+3*) 

s^r 


(2.39) 


It is seen from this equation that in the vicinity of 
resonance (w=oS^), the receptance is controlled by the first term 
in Eqn, (2.39) relating to the r-th mode. This term when drawn 
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in vicinity of a)= 0 }^ in Nyquist plot will be a circle like, 

Fig. 2, 3 (a) and second term (of series) giving the combined 
effects of all other modes is essentially constant through the narrow 
range of frequencies, 2, 3(b). Combined effect 

of two terms of the above equation is shown in Fig, 2.3(c), 
which is circle like . The greater the range of frequency over 
which the combined effect (off-resonant contribution) is constant, 
the greater will be the part of the curve coinciding with the 
circle. The circle is displaced by an amount determined by the 
combined effect of all other modes. 


When once the modal circle is fitted using mobilities 
in case of viscous damping and receptances in case of structural 
damping, the next step is to determine the natural frequency of 
that mode. The natural frequency, in case of viscously damped 
system, is located where the real component is maximum on the 
modal circle as shown in Fig, 2.4(a). Similarly the natural 
frequency, in case of structurally damped system, is located where 
the imaginary component on the modal circle is maximum. 


After locating the natural frequency, the next step 
is to estimate the damping of the mode under consideration. 

Consider viscously damped case. The mobility expression 
by considering Eqn. (2.36), can be written as follows. 








(2.40) 
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Considering Fig. 2.2(a), one can express the angle x as 

tan T= tan(©/2) = .... (2.41) 

wO'tO^+7 ' ) 

Taking ©^= 90° (below natural frequency) and ©^= 90° (above 
natural frequency), Fig. 2.4(a) 

0i^(i3' 5^+7' ) 

2 -2 

>] ^ 

W O' 5^+7 ') 
a^ r ^ 


Adding above equations, one gets 


+ 7 '= 

Using Eqn, (2.28), the above equation becomes 




2 C. 


Similarly considering structurally damped case, 
can be written as 


(2.42) 

Eqn. (2.38) 


ajk(^)= 


rg^k) . 

(w^-w^)^+Ow|+7)^ ^ 




(2.43) 
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Fig. 2.4 Propeties of modal circle. 
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Considering Fig. 2, 2(b), one gets the angle t as 


tan T = 


/3c3^+7 

r 

“-2 2 
w -w 
r 


or 


-2 2 
CO -<0 

tan(e/2) = tan(90-T) = — ^ 

3c3‘^+-r 

r 


(2.44) 


Taking 9j^=90° (below C^) and = 90° (above 5^) 
as shown in Fig. 2, 4(b) 


1 = 


05^4.7) 


1 = 


(352+7) 


Adding above equations, one gets 


2(352+7) = (o2-co2 


Using Eqn. (2.31), the above equation becomes 


n 

r 


252 

r 


« • • • ( 2 . 45 ) 


After determining natural frequency (5^) and damping 
(^^ or b^) at a given mode, the next step is to estimate the 
eigenvector associated with that mode. This is done by 

linking the diameter of modal circle with other two quantities 
viz, natural frequency and damping. 
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Consider Eqn. (2.4o) for viscously damped system 
substituting Oi=C^, one gets mobility as 




(w) 


C0=03 


r ,1 r k 


( 2 .^ 6 ) 


On observation, the above Eqn, contains only real 
part which is shown in Fig. 2, 5(a). The magnitude of 
co=w becomes the diameter (_D,, ) of the modal circle as the 

r JT jK 

Eqn. (2.40) doesn’t contain off -resonant terms. Rewriting above 
equation one gets 


r 


0 





(e-ct7-) = 2 5^ 




.. (2.47) 


Considering the point response at w=g 3^ one can 

compute (i.e. at u)=o> 2 , 5^ one can 

compute 2^1’ 3^1’ ^12’ ^1n^* 

first row of the normalized eigenvector matrix [ 0 ] is known. 

From transfer response y >|2 one can compute ^0^.^02 

(i.e. 0^^ 021 ^21 i^nown. Similarly coefficients 

^22’ ^^23’-**» ^2n computed from transfer response y^ 2 

w=S 2 , Continuing this, the full normalized 

eigenvector matrix [ 0 ] is determined. 

Simiarly the modal quantities can be related with the 
diameter of modal circle in case of structurally damped systems 


as follows : 
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Consider Eqn. (2. A3) for structurally damped case. 
Substituting Oii=w^ one gets, receptance as 




= -a 




05^+ T) 


(2.48) 


On observation, one can knov/ that the above equation 

contains only imaginary quantity, which is shown in Fig. 2.5(b). 

The magnitude of w=a)^ is equal to the diameter 

( D., ) of the modal circle as the Eqn. (2 A3) doesn’t contain 
r jK 

any terms pertaining to off -resonant contribution. Rewriting 
above equation, one gets 




= 

r^k 




(BCJ+t) » 


D., 

r ok r r 


co=w 


(2.49) 


Normalized eigenvector matrix ]is determined as 

explained earlier for viscous damping case. 

After knowing the modal matrix [0 ]» mass matrix [M ] 
and stiffness matrix [ K ] are obtained from Eqn. (2.22) as 


[M]=: [d [i ]-■' 

[K]- [d ]■’’ ["S^J 


... (2.50) 
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CHAPTER - III 


RESULTS AND DISCUSSION 


The first part of the present chapter discusses some 
of the assumptions made in the present analysis , In the latter 
part, the results obtained for the discrete and continuous 
systems with the structural /viscous damping are discussed. 

3.1 MODELLING METHODOLOGY 

In this section some of the simplifications made for the 
present analysis are explained below. 

The Eqn, (2.19) gives the relationship of the damping 
matrices with the mass and stiffness matrices. For the present 
analysis, the damping matrices are considered to be proportional 
to the stiffness matrices only. Hence, 

[C ]= ^'[k ] 

[H]=^[K] 


(3.1) 


With this assumption, 7=0 . The Eqns . (2.28) and 
(2.31) become 

.... (3.2) 


= 




As the experimental values of the modal response were 
not available, the receptance /mobility response information was 
generated from theoretical models (discrete and finite element 
models) . 
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3.2 RESULTS AND DISCUSSIONS 

Tables 3,1 -3. 1^- document the actual and identified 
values of mass and stiffness matrices for the cases studied. 

Figure 3.1 shows the two DOF system studied and 
Table 3,1 gives its actual mass and stiffness matrices. 
Theoretical receptance for structural damping ratio of 0 .01 
and 0 ,05 is calculated using Eqn. (2.38a), Mass and stiffness 
matrices are calculated from these using circle fit method. 

A comparision of results shows the close agreement between the 
actual and identified values . The identified mass is found to be 
almost diagonal, with the off-diagonal terms being small in 
comparision with the diagonal terms. The results for the case 
with a damping of 0.05 are closer than that of 0.01, as the 
response peaks are very well separated for the former case, thus 
making it easier to accurately select the range. This improves 
the quality of curve fit, thus making system identification 
more accurate. 

Figure 3.2 shows an axial rod fixed at one end. First it 
was modelled as two finite elements and then as three finite 
elements. Consistent mass matrices were used. These matrices are 
given in Table 3. 2-3. 3. Structural damping factor taken were 0.01 
and 0.05. Theoretical receptance was calculated using Eqn. (2.38a). 
Mass and stiffness coefficients were calculated from this using 
circle fit method. Identified values are given in Table 3. 2-3. 3. The 
identified values agree well with the actual values . 
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The beam with one end fixed is shown in Fig. 3,3. Here 
also it was modelled first as two finite elements and then as 
three finite elements with 2 DOF per node. Consistent mass matrices 
were used. Theoretical receptance was generated using Eqn. (2.38a). 
The mass and stiffness values were identified using circle fit 
method. Both assumed and identified values are given in 
Table 3. 4-3. 7 for structural damping factor 0.01 and 0.05. 

It is observed that the identified values agree closely with 
actual values . 

Figure 3.4 shows L-Frame structure fixed at node 3. It 
was modelled as two elements with three DOF per node. Consistent 
mass matrices were used and theoretical receptance was generated 
using Eqn. (2,38a). Table 3. 8-3. 9 show the identified and 
actual values for structural damping factor 0,01 and 0,05, It is 
observed that the identified values are closer to actual values. 

Figure 3.5 shows F-Frame structure fixed at node 5. 

It was modelled as 4 elements with three DOF per node. Consistent 
mass matrices were used. Theoretical receptance was generated 
using Eqn. (2.38a). The identified mass and stiffness values are 
shown in Table 3.10. The identified values are quite closer 
to actual values . 

Some of the cases discussed earlier, were also analysed 
for proportional viscous damping. Consistent mass matrices were 
used. Theoretical receptance was generated using Eqn. (2.38a)* 



27 


Later on, mobility was obtained from receptance for analysis . 

The modal parameters were extracted using circle fit method. 

The various cases analysed are discussed below. 

Table 3,11 shows the actual and identified values 
for a two DOF discrete system, Fig. 3.1 , The identified 
values agree well with actual values . 

Table 3.12 shov\/s the actual and identified values for 
an axial rod with two elements. Fig. 3.2(a). The identified 
values agree well with actual values . 

Table 3,13 shows the identified and actual values 
for a beam with two elements. Fig. 3.3(a), The identified values 
agree well with actual values . 

Table 3.14 shows the identified and actual values for 
L-Frame, Fig. 3.4. The identified value agree well with 
actual values . 

3.3 CONCLUSION 

From the present analysis it can be concluded that 
circle fit method gives fairly accurate identification of modal 
parameters for discrete as well as continuous structures with 
viscous /structural proportional damping. 



in'] = 1.0 kg ; k-] = 15.0 kN/m 
m 2 = 2.0 kg 5 k 2 = 20.0 kN/m 


Fig. 3.1 Two DOF discrete system 
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Table 3.1 Comparision for a discrete system (2x2) 

(structural damping) 


n=0.01 


Actual 

values 

Identified values 


1 .0000 

0 .0000 

1 .0148 

-0 .0009 

Mass 





matrix 

0 .0000 

2.0000 

-0 .0009 

2.0496 

Stiffness 

3.5000 

-2 .0000 

3.5464 

-2.0496 

matrix 





X 1(7^ 

-2,0000 

2 .0000 

-2.0496 

2.0642 


Ti= 0.05 



Actual 

values 

Identified values 

Mass 

1 .0000 

0 .0000 

0.9992 

-0.0039 

matrix 

0 .0000 

2.0000 

-0.0039 

2.0132 


3.5000 

-2.0000 

3.5079 

-2.0061 

Stiffness 





matrix 

. "^4 

X 10 

-2 .0000 

2 .0000 

-2.0061 

2.0074 


(a) Two elements 


L 600 mm 

A 625 mm^ 

p 7800 kg/mm 

E 207 GPa 



12 3 4 

(b) Three elements 


Fig. 3.2 Axial rod, finite elements. 
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Table 3.2 Comparision for an axial rod(2x2) (structural damping) 


n = 0.01 


Actual 

values 

Identified values 

Mass 

5.1071 

2.5535 

5.0923 

2.5305 

matrix 

2.5535 

1 .0214 

2.5305 

1 .01 81 

Stiffness 

4.1250 

- 4.1250 

4.1231 

- 4.1299 

matrix 

X ia® 

- 4.1250 

8.2499 

- 4.1299 

8.2455 


T ) = 0.05 


Actual 

values 

Identified values 


5.1071 

2.5535 

5 .1066 

2.5365 

Mass 





matrix 

2.5535 

1 .0214 

2.5365 

1 . 01 33 


4.1250 

- 4.1250 

4.1172 

- 4.1164 

Stiffness 





matrix 

X 10® 

- 4.1250 

8.2499 

- 4.1164 

8 . 21 85 


Table 3 .3 .Comparision for an axial rod (3x3) (structural damping) 
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1 2 3 

(a) Two elements 


L 600 mm 
A 625 mm^ 

I 32552 mm^ 
p 7800 kg/m^ 
E 207 GPa 



12 3 4 

(b) Three elements 


Fig. 3.3 Beam finite elements. 



Table 3,4 Comparision for a beam (4x4) (structural damping) 
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Table 3.5 Comparision for a beam (Ax4) (structural damping 
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Table 3.6a Comparision for a beam (6x6) (structural damping 
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contd . 



Table 3.6b Comparision for a beam (6x6) (Structural damping 
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1 

1 


77777: 



L-i 3 00 mm 

L 2 300 mm 

A 625 mm^ 

I 32552 mm^ 

p 7806 kg/m^ 

E 207 GPa 


Fig. 3.4 L- Frame structure . 



Table 3.8a Coniparlsion for L-Frame (6x6) (structural damp! 
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Table 3,Qb Comparision for L-Frame (6x6) (Structural damping 


42 




o 

in 

m 

o 

in 

00 


T“ 

4 

4 

m 

CM 

CO 



o 


o 

in 

4* 

r- 


o 

4 

O 

4 

m 

T- 



o 

o 

o 

o 

o 

O 


o 

O 

o 

O 

o 

o 



o 

o 

o 

o 

o 

O 


o 

O 

o 

o 

o 

o 



« 

* 

• 

* 

* 

• 


« 

% 

• 

• 

• 

• 



o 

o 

o 

9 

9 

o 


9 

o 

o 

9 

9 

o 



o 

o 

tn 

o 

in 

in 



in 

o 

n- 

M3 

CM 



o 

o 

<t 

o 

CM 

4 

00 

n- 

CM 

CM 


m 


03 

o 

m 

o 

o 

4” 

O 

*o 

o 


o 

O 

rn 

o 


o 

o 

o 

o 

m 

O 

T“ 

o 

o 

o 

O 

o 

o 


KD 

• 

* 

• 

• 

« 

• 

X 

* 

• 

• 

• 

• 

• 


T” 

o 

9 

9 

o 

4- 

9 

X 

9 

9 

9 

o 

4 

9 


X 







p 























X 







p 








•H 







cd 








u 







s 








•p 















cd 







CO 





r- 



s 

m 

O 

o 

in 

O 

O 

CO 

o 

M3 

T“ 

m 

m 



CM 

O 

o 

CM 

O 

in 

(U 

o 

O 

o 

o\ 

CM 

4 


m 

T~ 

o 

o 


O 

o 


o 

o 

o 

CM 

o 

O 


to 

ro 

o 

o 

m 

O 

o 

p 

m 

o 

o 

m 

o 

o 

o 

0) 

« 

« 

# 

• 

• 

« 

p 

♦ 

• 

• 

« 

• 

• 

o 

p 

i 

o 

o 

•4- 

o 

9 

•H 

P 

4" 

1 

o 

o 

4 

o 

9 


ch 







to 







II 

•H 














P 







Tj 







p 

W j 







0) 















I 








r— 1 







P 1 








Cd ! 







*ri 















i 








p ■ 







^ \ 








o 







(D ' 







1 

C ; 







P ' 















M j 






4 



o 

in 

cn 

O 

in 

in 


T" 

tn 

cn 

r- 

Q 

1 


o 


o 

o 

4” 

o 


o 

4 

o 

O 

CM 

O 



o 

o 

o 

o 

o 

o 


o 

O 

o 

o 

O 

o 



O 

• 

o 

• 

o 

• 

o 

« 

o 

• 



o 

• 

o 

• 

o 

• 

o 

• 

O 

* 

o 

« 



o 

o 

o 

o 

9 

o 


9 

o 

o 

o 

9 

o 


■ 

o 

o 

m 

o 

o 

m 


in 


m 

MO 

in 

4 



o 

o 


o 

o 

4 

! 

o 


4 

o 


4 



o 

m 

o 

o 

m 

o 


o 

CM 

O 

o 


O 



o 

• 

q 

o 

• 

o 

« 

o 

• 

o 

♦ 


o 

• 

O 

• 

o 

• 

o 

• 

o 

« 

o 

* 



o 

o 

o 

o 

9 

o 


9 

o 

o 

o 

9 

o 



in 

CM 

o 

o 

o 

o 

in 

CM 

O 

o 

o 

o 


m 

o 

in 

o 

r- 

O 

o 

o 


o 



T- 

m 

o 

o 

• 

o 

o 

o 

• 

o 

m 

• 

■4- 

1 

o 

o 

• 

o 

o 

o 

• 

o 


o 

m 

• 

4 

,o 

o 

9 

o 

o 

9 

0 
m 

• 

4 

1 

O 

o 

* 

9 

o 

o 

9 
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Table 3.10a Comparision for F-Frame (12x12) (Structural Damping 
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Table 3.11 Comparision for a discrete system (2x2) 

(viscous damping) 


iB ' = 1 .0 X l 5 ^ 


Actual 

values 

Identified values 


1 .0000 

0.0000 

1 .0303 0.0473 

Mass 




matrix 

0 .0000 

2.0000 

0.0473 1.9397 



3.5000 

- 2 .0000 

3.5070 

-1 .9755 

Stiffness 





matrix 





X 10 ”-^ 

- 2.0000 

2.0000 

-1 .9755 

1 .9704 
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Table 3,12 

Comparision for an axial rod 
(viscous damping) 

(2x2) 



j3'= 

1 .0 X 10^ 




Actual values 

Identified values 


5.1071 

2.5535 

5.0786 

2.5252 

Mass 





matrix 

2.5535 

10.21 43 

2.5252 

10.1591 


4.1 250 

-4.1250 

4.1122 

-4.1188 


Stiffness 

matrix 

X 10-S 


-4.1250 


8.2500 


-4.1188 


8.2246 


Table 3.13 Comparision for a beam (4x4) (viscous damping 
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Table 3 Comparision for L-Frame (6x6) (viscous damping) 
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Table 3.14b Comparision for L-Frarae (6x6) (viscous damping) 


52 




o 

LTs 

in 

o 

m 

CO 


CA 

r- 

A 

A 

r- 

00 



o 

-4“ 

o 

m 


r~- 


O 

<f 

o 

A 

o 

T™ 



o 

O 

o 

o 

o 

O 


O 

o 

o 

O 

o 

o 



o 

O 

o 

o 

o 

O 


o 

o 

o 

O 

o 

o 



• 

« 

• 

* 

• 

« 



« 

* 

* 

• 

« 



o 

o 

o 

9 

9 

o 


o 

o 

o 

9 

o 

o 









00 















kD 























CO 

o 

o 

m 

O 

tn 

m 

X 

<r 

A 

MO 

A 

CO 

r- 


*o 

o 

o 

<r 

o 

CM 



T" 

X“ 

A 

A 

00 

o 


V" 

o 

ro 

o 

o 

<r 

o 

X 

r“ 

— 

O 

T— 

o 

o 


X 

o 

O 

o 

o 

fA 

o 

•H 


O 

o 


CM 

o 



« 

• 

• 

• 

• 

• 

u 

« 

* 

• 

• 

• 

• 


X 

•H 

o 

9 

9 

o 


9 

•p 

cd 

O 

o 

o 

9 

A 

o 


u 







a 






























CO 








B 







to 















(1) 








w 







a 







»o 

CO 














Y- 

0) 















a 

m 

O 

O 

m 

o 

o 


CA 

n- 

A 

CO 

A 

A 

X 


C\J 

o 

o 

CM 

o 

lA 

4J> 


r- 

T" 

CM 

A 

A 




o 

o 

-4- 

o 

O 

CO 

r- 

o 

o 

Q 

r“ 

O 

o 

•H 

r<^ 

o 

o 


o 

o 


A 

o 

o 


IT" 

O 

• 

•P 


• 

• 

* 

• 

• 

xi 

♦ 

« 

• 

• 

• 

• 


CO 

-4* 

1 

o 

o 


o 

9 


-4* 

1 

9 

9 


9 

9 

11 

rH 







Cp 








cd 







•H 







ca 

:z5 







-P 








*p 







d 








o 







0) 








< 






















HH 









o 


cn 

a 

cr\ 

m 


T" 

A 

A 

A 

MO 

A 



1 o 

-4- 

o 

o 

<r 

o 




O 

r** 

A 

O 



o 

O 

o 

o 

o 

o 


o 

O 

O 

O 

O 

O 



o 

o 

o 

o 

o 

o 


o 

O 

O 

O 

O 

O 



• 

* 

• 

• 

« 

% 


• 

• 

• 

• 

• 

• 



o 

o 

o 

o 

9 

o 

i 

1 

o 

o 

o 

9 

o 

o 



o 

o 

in 

o 

O 

A 


1 


A 

r- 

A 




o 

o 


o 

O 

-4" 



o 


A 

r" 




o 


o 

o 

m 

O 

j 

o 

A 

O 

o 


o 



o 

o 

o 

o 

o 

O 


o 

o 

o 

o 

O 

o 



« 

• 

• 

• 

• 

• 


• 

• 

• 

• 

« 

* 



o 

o 

o 

o 

9 

o 


o 

o 

o 

9 

o 

o 



in 

o 

o 

ir\ 

O 

o 


CA 

Y“ 

<r 

A 

<r 

A 



CM 

o 

o 

CM 

o 

o 


O 

I> 

r- 

-r* 

r- 

O 



^T" 

o 

o 

-T” 

o 

o 


r- 

o 

O 

r- 

T* 

O 



rr\ 

o 

o 


o 

o 


A 

o 

O 

A 

T“ 

O 



• 

• 

« 

• 

• 

41 


• 

• 

• 

« 

• 

« 



<r 

o 

o 

<r 

o 

o 


-4“ 

o 

o 


o 

o 






1 







t 





REFERENCES 


Kundra, T.K., 'Studies in Identification and Modification 
of Dynamic Mechanical Systems' , Ph.D. Thesis, Dept, of Mech. 
Engg., Delhi, April 1986. 

Vepa, K., 'Optimal Identification of Vibrating Structures', 
International Modal Analysis Conference and Exhibit (IMAC2), 
Orlando, Florida, Feb. 1984. - 

Ewins, D.J., 'Modal Testing : Theory and Practice', Research 
Studies Press Ltd., John Wiley and Sons Inc .#1984. 

Kennedy, C.C., and Pancu, C.D.P., 'Use of Vectors in 
Vibration Measurement and Analysis', J. Aeronautical Sciences, 
Vol. 14, No. 11, 1947, pp. 603-625. 

Bishop, R.E.D., and Gladwell, G.M.L., 'An Investigation into 
the Theory of Resonance Testing', Proc Roy Soc Phil Trans 
255(A) 241, 1963. 

Salter, J.P., 'steady State Vibration', Kenneth Mason Press, 
1969. 

Allemang, R.J., 'Experimental Modal Analysis Bibliography’, 
Proc. IMAC2, Feb. 1984. 

Mitchell, L.D., and Mitchell, L.D., 'Modal Analysis 
Bibliography - An update- 1980-1983 Proc. IMAC2, Feb. 1984. 



54 


9. Flannelly, W,G., McGarvey, J.H., and Berman, A., ' A Theory 
of Identification of the Parameters in the Equations of 
Motions of a Structure through Dynamic Testing Papers 
Presented at the Symposium on Structural Dynamics, University 
of Technology, Loughborough, U.K., March 23-25, 1970. 

10. Berman, A., ' Determining Structural Parameters from Dynamic 
Testing', The Shock and Vibration Digest, 7(1), Jan. 1975. 

11. Zaveri, 'Modal Analysis of Large Structures - Multiple 
Exciter Systems’ , Bruel and Kjaer, April 1985. 

12. Meirovitch, L., 'Elements of Vibration Analysis' , McGraw- 
Hill International Editions, 1986. 



